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Abstract 

We study a functional method to extract the V—A condensate of dimension 6 from 
a comparison of r-decay data with the asymptotic space-like QCD prediction. Our 
result is in agreement within errors with that from conventional analyses based 
on finite energy sum rules. 



* Supported by PROCOPE, EGIDE 
tEmail: sebu@lpm.univ-montp2.fr 

* Email: hspiesb@ thep . physik . uni-mainz .de 



1 Introduction 



Although QCD has been with us for three decades, the knowledge of the values of the 
various fundamental or effective parameters of the theory (with the possible exception 
of the coupling constant) such as quark masses and condensates is still astonishingly 
limited. The precise data on r-decay obtained by the ALEPH jT] and OPAL [2] col- 
laborations at CERN have offered an opportunity for new studies, which range from 
an extraction of the strange quark mass [Sj to the determination of various condensate 
parameters. Of particular interest was the extraction of the dimension-6 condensate 
[H El E3 El IE] which, in the chiral limit, determines, e.g., the K — > tttt matrix elements 
of the relevant electroweak penguin operators. 

It should be kept in mind, however, that the extraction of the condensate parameters 
of QCD constitute a so-called ill-posed inverse problem which is basically unstable 
with respect to errors in the input data. An example is the popular method used to 
obtain condensates from experimental spectral functions by QCD sum rules (finite- 
energy sum rules, FESR). This procedure corresponds to an analytic continuation from 
a finite contour (the part of the positive real axis on which the data are given) which is 
known to be notoriously unstable with respect to data errors. The extraction of QCD 
condensates requires therefore a carefully chosen stabilization mechanism. In the case 
of QCD sum rules, this is achieved by an implicit and rather ad hoc assumption that the 
series of the operator product expansion essentially breaks off after a finite number of 
terms. Amazingly, a careful analysis of finite energy sum rules [7j of the chiral spectral 
function shows remarkable stability with respect to variations of the duality interval. 
To investigate further the reliability of the extraction of QCD parameters via sum rules, 
it would be prudent to develop alternative methods. 

In this letter we study a functional method [31101 which allows us to extract without 
prejudice the condensates from a comparison of the time-like data with the asymptotic 
space-like QCD results. We will see that the price to be paid for the increased credibility 
are possibly larger errors in the values of the extracted parameters. 

2 QCD condensates 

We consider the polarization operator of hadronic vector and axial-vector charged cur- 
rents, J M = Vfj, = uj^d and J^ = A (i = wy^d, 



The conservation of the vector current implies Tl v = 0. 

The spectral functions are related to the absorptive part of the correlators 
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and can be measured in hadronic r-decays. We consider specifically the V — A compo- 
nent which is related to the branching ratios of r decays through 

B{j — > v T + hadrons, V — A) 
R T y-A = ^7 ; — t (3) 
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Here, V u d is the weak mixing CKM-matrix element, \V u d\ = 0.9752 ± 0.0007, the r 
mass is denoted m T = 1.777 GeV and Sew = 1.0194 ± 0.0040 accounts for electroweak 
radiative corrections ^T]. The spin-0 axial vector contribution ao(s) is dominated by 
the one-pion state, ao(s) = 2ir 2 f%5(s — m%), with the 7r-decay constant f n = 0.1307 
GeV. Its contribution in the last term of Q is tiny 

AR rV _ A \ ~ 24tt 2 ^-^ ~ 0.0074 (4) 

and will be neglected. The contribution of the pion pole to the first term in Eq. © 
is well identified in the data and concentrated at low s ~ m\] thus it can be removed 
from the data and taken into account explicitly without introducing sizable additional 
uncertainties. The experimental data PJE]* are given by binned and normalized event 
numbers related to the differential distribution dR T y_A/ds and can therefore be viewed 
clS db measurement of the function 

uj V -a(s) = v^s) - a^s) - a (s). (5) 

The (V — A) correlator is special since it vanishes identically in the chiral limit 
(m q = 0) to all orders in QCD perturbation theory. Renormalon ambiguities are thus 
avoided. Non-perturbative terms can be calculated for large \s\ by making use of the 
operator product expansion (OPE) of QCD 



n ») = Ei5( 1+c v) < 6 > 



D>A 



where O v D are vacuum matrix elements of local operators of dimension D (so-called 
condensates). Their contribution is known up to dimension 8 and read, at leading order, 



0\~ A = (m u + m d )(qq) = -f*ml, (7) 
Ol~ A = - 3 -^a s (qq)\ (8) 
O v & " A = Aixa s i{qG aP G^q). (9) 



The last two results hold in the vacuum dominance approximation. The numerical 



value of 0\ A is very small and this condensate can be neglected in our analysis. 



*We use the ALEPH data because of their smaller experimental errors. 
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The next-to-leading-order corrections to 0$~ A have been calculated in [TJl[T3j. They 
depend on the regularization scheme implying that the value of the condensate itself is 
a scheme- dependent quantity. Explicitly, 
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anticommuting 75 [T5] . 
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The renormalization scale fi 2 is conveniently chosen to be —s. 

The result for Og~ A , Eq. (|9*|). is taken from ^6J. It involves a quark-gluon condensate 
for which various estimates exist. The typical scales determining the condensates are 
around 300 MeV, e.g. (qq) ~ (250 MeV) 3 , (a fl /7r)(G 2 ) ~ (300 MeV) 3 . Assuming a 
similar scale for the condensate entering we expect 0%~ A to be of order 10 -3 

GeV 8 . This is small enough so that the OPE makes sense. If 0\~ A would be lar ger, 
radiative corrections to higher-dimension condensates would mix significantly with the 
lower-dimension condensates through their imaginary parts. There exist a number of 
QCD sum rule extractions of the value of the D = 8 condensate. They range from 
(-7-5±4:o) • 10- 3 GeV 8 to (4.4 ± 1.2) • 10~ 3 GeV 8 @|. A recent conservative estimate 
[7j is 8 = (-1.0±6.0) • 10~ 3 GeV 8 . This value corresponds to a scale of about 400 MeV 
which is comparable to Aqcd- The variation of these results represents the ambiguities 
inherent in the QCD sum rule approach. In the next section we shall present our 
alternative rigorous functional method which allows us to extract the condensates from 
a comparison of the data with the asymptotic space-like QCD results in an unambiguous 
way. 



3 An L 2 norm approach 

We consider a set of functions F(s) (where F(s) relates to Hy^(s)) expressed in terms 
of some squared energy variable s which are admissible as a representation of the true 
amplitude if 

i) F(s) is a real analytic function in the complex s-plane cut along the time-like 
interval Tr = [s ,oo). The value of the threshold so depends on the specific 
physical application (s = (2m jr ) 2 for H v , s = m 2 , for Ha)- 

ii) The asymptotic behavior of F(s) is restricted by fixing the number of subtractions 
in the dispersion relation between F(s) and its imaginary part along the cut 
f(s) = ImF(s + i0)\ S £r R (for Hyt A \s) no subtractions are needed): 

F(s) = - !°° ^-dx. (12) 

7T /.„ X — S 
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We have two sources of information which will be used to determine F(s) and f(s). 
First, there are experimental data in a time-like interval r exp = [s , s max ] with sq > 
for the imaginary part of the amplitude. Although these data are given on a sequence 
of adjacent bins, we describe them by a function f exp (s). We assert that f exp is a real, 
not necessarily continuous function. The experimental precision of the data is described 
by a covariance matrix V(s, s'). 

On the other hand, we have a theoretical model, in fact QCD. From perturbative 
QCD we can obtain a prediction for the amplitude in a space-like interval Fi = [s2, si]. 
This model amplitude Fq CD (s) is a continuous function of real type, but does not 
necessarily conform to the analyticity property i). Since perturbative QCD is expected 
to be reliable for large energies, we expect that there is also useful information about 
the imaginary part of the amplitude provided that \s\ is large, i.e. we can also use 
/qcd(s) = ImFQ CD (s + z'0)| sg ( SmaxiOO ). In order to compare the true amplitude with 
theory, we can therefore split the integral in the dispersion relation ([12)1. 

Fw-if m dx= ir-m dXt (13) 

7T /, X — S 7T X — S 

and test the hypothesis whether the left-hand side can be described by QCD. 

We also need an a-priori estimate of the accuracy of the QCD predictions. This will 
be described by a continuous, strictly positive function cxl(s) for s G Tl which should 
describe errors due to the truncation of the perturbative series and the operator product 
expansion and is expected to decrease as \s\ — > oo and diverge for s — > 0. In the case 
of IV-A which does not have perturbative contributions, we will take the contribution 
of the dimension D = 8 operator as an error and use <Jl{s) = 0$/s 4 with 0$ in the 
order of 1CT 3 GeV 8 . If the perturbative part dominates, as is the case for the individual 
vector or axial vector correlators, the last known term of the perturbation series could 
be used as a sensible estimate of the error corridor. 

The goal is to check whether there exists any function F(s) with the above analyt- 
icity properties, the true amplitude, which is in accord with both the data in r exp and 
the QCD model in T^. In order to quantify the agreement we will define functionals 
Xtif] an d XrU'] using an L 2 norm. For the time-like interval we simply compare the 
true amplitude f(s) with the data and use the covariance matrix of the experimental 
data as a weight function: 

dx / dx' V~\x, x') (f(x) - f exp (x)) (f(x') - / exp (x')) • (14) 

o J So 

Experimental data correspond to cross sections measured in bins of s, so that we can 
calculate this integral in terms of a sum over data points. The ALEPH data which we 
use are given for 65 equal-sized bins of width As = 0.05 GeV 2 between and 3.25 GeV 2 . 
x\ given in ([14)1 is in fact the conventional definition of a \ 2 ail d has a probabilistic 
interpretation: for uncorrelated data obeying a Gaussian distribution we would expect 

t\Ve do not exclude the case S2 — > — oo. 
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to obtain x\ — N, where N is the number of data points. Since experimental data at 
different energies are correlated, we instead expect 



In order to define a measure for the agreement of the true function f(s) with theory, 
we use the left-hand side of (fT3|) which is well-defined and expected to be a reliable 
prediction of QCD in the space-like interval for not too small \s\. This expression can 
be compared with the corresponding integral over the true function. Thus we define 

x L \f] = N J t wl{x) [f qcd (x) --j dx - -J -^rr^ dx ) dx 

(16) 

where wl is the weight function for the space-like interval and identified with 1/<j 2 l (s). 
The integral is normalized to unity for the case where the difference within parentheses 
saturates the error a L . 

In order to find the true function f(s), we can combine the information contained 
in Xr-i flSj) , and \\i (HHJ) in the following way (HIHD!- We fix 

X R if} = xl<xl P , (17) 

and minimize xi'- 

Xl[f}^ least (=xl M n)- (18) 
These conditions are equivalent to finding the unrestricted minimum of the functional 

m=xi[f]+nx 2 R [f] 

where \i is the Lagrange multiplier, which will be found later. The solution of the 
condition x\\f \ — Xo w ^ b e denoted by f(x;fM): 



max 



dx / dx'V 1 (x, x') [f(x; //) - f cxp (x)] [f(x'; //) - f cxp (x')] = xl ■ 
To this end we require the Frechet derivative of T to be zero 

oF[f, Y] = lim n 5 = , 

a^o da 

for any function Y. This leads to the following integral equation for the imaginary part 

= /«p(x) + - / dx'V(x,x') dx"w L (x")F QC T)(x")— 

" 'so Jv ' " - " 
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(x' - y)(x" - y) 

+ A / dx'JC(x,x')f(x';n), (19) 
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where A = l//i and 



K(x, x) 




dx"V(x, x") \ dy 



wl{v) 



(x' — y)(x' 
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Eq. (fTHI) is a Fredholm equation of the second type which is stable against variations of 
its input. At this stage we should notice that if one had claimed that in the space-like 
region the function F(s) was given by some analytic expression (e.g., by some few QCD 
terms), this would be equivalent to saying that x\ vanished identically. But then, from 
the definition of the functional JF and the vanishing of its Frechet derivative, it follows 
that fi = 1/A is zero which will turn the integral equation (|T9|) into a Fredholm equation 
of the first kind which is known to be unstable. 

The integral equation will be solved numerically by expanding f(s) in terms of Leg- 
endre polynomials. The algorithm to determine an acceptable value for the condensate 
is then the following: 

i) For a fixed value of xl — xLp we determine the solution (fT^j) and calculate 
the corresponding value of xVlf] as a function of the condensate a s {qq) 2 . The 
Lagrange multiplier /i is determined by iteration such that the condition xWf] = 
Xo is fulfilled. 

ii) We minimize this Xzlf] with respect to a s (qq) 2 and call the minimal value ximm 
and the corresponding a s (qq) 2 is the value for the condensate we are looking for. 

Hi) We determine the error on a s (qq) 2 by solving x|( a s(^) 2 ) — xi m in + 1- 

4 Numerical results and discussion 

A typical situation resulting from this algorithm is shown in Fig. ^ The left part of this 
figure shows x\ which has the expected quadratic dependence of a s (qq) 2 . The values of 
a s (qq) 2 corresponding to the minimum of x\ are listed in the Table below for various 
choices of Os as discussed above. The regularized function shown in the right part of 
Fig.Hfollows nicely the data points, except at large s. Here the experimental errors are 
large and hence, as it should happen, the regularizing effect by means of the functional 
(fTlj) is not as effective. 

N f=3 

In the numerical evaluation we have used the NLO expression for a s (s) with = 
0.326 GeV. The result for a s (qq} 2 is not sensitive to changing A^ -3 within the present 
experimental error ±0.030 GeV. For the evaluation of x\ we have restricted the range 
of integration within limits S2<s<si<0. We checked that our result is insensitive 
to changes of s 2 as soon as its absolute value is chosen larger than (9(100) GeV 2 . Since 
the error channel defined by Os/s diverges for s — > 0, one could, in principle, choose 
the upper limit si = 0. Numerical instabilities require a non-zero value. We observe 
a well-defined plateau for the result for a s (qq) 2 as a function of si between —1.0 and 
—0.5 GeV 2 and quote the values for si = —0.7 GeV 2 . 
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function 
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: A typical result for xt as a function of u s (qq) 2 (left) and the regularized 
compared with data pQ (right). We have chosen Os = 1CT 3 GeV 8 and cq = 



The values for a s (qq) 2 given in the table translate into values for the condensate 
0\~ A according to Eq. (fTUj) . In order to compare with other results from the literature 
we use a s (s) = 0.6 at the scale s — 1 GeV 2 . For 8 = 1.0 x 10 -3 GeV 8 we obtain 



V-A 



(-0.0020 ± 0.0014) GeV 6 for c 6 = § , 
(-0.0015 ± 0.0009) GeV 6 for c Q = ^§. 

These results can be compared with the lowest-order vacuum saturation expression 

32vr 
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-0.0013 GeV 3 



where we used (qq) = —0.014 GeV 3 . On the other hand, analyses based on finite energy 







a s {qq) 2 for c 6 = § 


a s (qq) 2 for c 6 = 2 § 


8 = 1.0 x 10~ 3 


GeV 8 


1.6 ± 1.0 


1.1 ±0.6 


1.25 x 10~ 3 


GeV 8 


1.6 ±1.1 


1.1 ±0.8 


1.5 x 10~ 3 


GeV 8 


1.6 ±1.2 


1.1 ±0.9 



Table 1: Results of the determination of a s (qq) 2 (in units of 10~ 4 GeV 6 ) for the two 
choices of Cq in (jllj) and with different values for Og to fix the error channel in the 
space-like interval. 
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sum rules [H 03 El IE] typically find results 

6l~ A = (-0.004 ± 0.001) GeV 6 , 

which are not inconsistent with our number. The fact that we find agreement within 
errors is not trivial. Since our approach is based on less assumptions we may conclude 
that the sum rule results with their relatively small errors are indeed trustworthy. 
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